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Abstract 

We use numerical and asymptotic techniques to study the stability of a two- 
phase air/water flow above a flat porous plate. This flow is a model of the boundary 
layer which forms on a yawed cylinder and can be used as a useful approximation 
to the air flow over swept wings during heavy rainfall. We show that the interface 
between the water and air layers can significantly destabilize the flow, leading to 
traveling wave disturbances which move along the attachment line. This instability 
occurs for lower Reynolds numbers than is the case in the absence of a water layer. 

We also investigate the instability of inviscid stationary modes. We calculate 
the effective wavenumber and orientation of the stationary disturbance when the 
fluids have identical physical properties. Using perturbation methods we obtain 
corrections due to a small stratification in viscosity, thus quantifying the interfacial 
effects. Our analytical results are in agreement with the numerical solution which 
we obtain for arbitrary fluid properties. 
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puter Applications in Science and Engineering (ICASE), M/S 132c, NASA Langley Research Center, 
Hampton, VA, 23681-0001. This work was also supported by the Science and Engineering Research 
Council. 
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1 Introduction 


The laminar flow over an infinitely long cylinder can become unstable as the Reynolds 
number increases. When the axis of the cylinder is inclined at an angle relative to 
the free stream, the developed three-dimensional mean flow can be separated into two 
components, one lying in a plane normal to the axis, the other parallel to the generators 
of the cylinder. Small amplitude disturbances to the flow can take the form of Tollmien- 
Schlichting waves, crossflow vortices, or Taylor Gortler vortices (if there are regions of 
concave curvature). 

The flow over a swept cylinder has been studied in detail, primarily because of its 
important application and relevance to the boundary layer which forms on the surface of 
swept-wing. Understanding the mechanisms of flow instability for this model, can lead 
to significant development of methods used in the reduction of laminar to turbulent flow 
transition. 

The model we use in this paper, is a classical Ileimenz stagnation point flow, together 
with a superposed non-zero component of velocity parallel to the axis. The equations 
governing the flow are written in cartesian coordinates (see Prandtl [24]). The velocity 
component parallel to the axis of the cylinder can be determined independently by 
decoupling the momentum equations. The relevance of this solution to the realistic flow 
which forms on a swept wing is discussed in section 3.2. 

Using linear stability theory, Hall, Malik & Poll [12] calculated critical Reynolds 
numbers for an infinite swept attachment line boundary layer. They examined the effects 
of both suction and blowing at the boundary. Surface suction can be used as an effective 
laminar flow control since it thins the viscous boundary layer and leads to a reduction 
in the local Reynolds number. In addition, the vorticity distribution is modified so that 
a more stable flow is established. Hall et al. obtained numerical and asymptotic results 
which clearly illustrate that even a small amount of suction can significantly stabilize 
the flow. Their results are in excellent agreement with the experimental investigations 
of Gaster [8], Pfenninger & Bacon [21] and Poll [22], [23]. These authors investigated the 
stability of the attachment lines on swept wings and swept cylinders to small disturbances 
of naturally occurring frequences. 

In 1986, Hall & Malik [11] extended their linear stability results to include the non- 
linear regime. The weakly nonlinear stability of this flow was examined using a Stuart- 
Watson expansion procedure. The primary motivation was to explain why experimental 
observations all correspond to modes near the lower branch of the neutral curve. Hall 
& Malik showed that apart from a small region near the critical Reynolds number, fi- 
nite amplitude solutions bifurcate from the upper branch when the Reynolds number 
is below the neutrally stable value elucidated from a linear stability analysis (subcrit- 
ical). Equilibrium states associated with the upper branch are not therefore observed 
experimentally, since these solutions are unstable. 

In addition these authors used numerical methods to integrate the time-dependent 
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Navier-stokes equations which govern the fully nonlinear problem. Using a Fourier- 
Chebyshev spectral method Ilall k Malik found the existence of supercritical finite- 
amplitude states near the lower branch of the neutral curve. 

Recently there has been much interest in the aerodynamic penalties associated with 
adverse weather conditions on aircraft flight. In a review of recent studies into the effects 
of heavy rain during take-off and landing, Dunham, Dunham k Bezos [6] showed that 
short duration, heavy precipitation can result in a premature loss of lift of 15 — 20% and 

an increase in drag coefficient of up to 20%. 

The exact mechanisms which cause these significant flight characteristics are not 
clearly understood. One possible explanation is that the presence of a thin water layer 
on the wing surface leads to a loss in stability of the laminar air flow. The growth 
of small disturbances in either the water or air layers could then promote transition. 
In this paper, we consider the interfacial stability of an infinite swept attachment line 
boundary layer consisting of a water layer and air flow above. The interface which 
separates these two viscous fluids may be susceptible to instability of the form first 
discovered by Yih [28]. Yih showed that long wavelength disturbances to plane Couette 
and plane Poiseuille flow of two immiscible liquids of different viscosities and densities 
can lead to a mode of instability which is a due entirely to the discontinuity in the fluid 
properties. The growth rate of the interfacial deflection approaches zero asymptotically 
as the viscosities of the two fluids become equal, hence this mode is not operational for 
a single fluid. The relative depths of the fluid layers is a crucial factor in characterizing 
the flow instability. 

Since Yih's work, there have been numerous investigations of interfacial instability 
which have important applications in many situations. For example, Blennerhassett [2] 
showed that the interfacial instability of air flow over water can lead to the generation of 
finite amplitude waves. The effects of surface tension and gravity have been quantified 
in a variety of numerical and analytical studies which consider short, moderate and long 
wavelength perturbations to the basic state (see Hooper k Boyd [15], [16], Hooper [14] 
and Renardy [25]). 

In this paper we quantify the effect of interfacial stability on the air flow over a 
swept wing during heavy rainfall. Using the model described above we obtain an exact 
solution of the Navier-Stokes equations which govern the viscous two-phase flow. The 
domain consists of two separate regions. In the upper region of the boundary layer we 
have a two-dimensional stagnation point air flow together with a superposed crossflow 
component (due to the angle of inclination to the free-stream). Below the air, is a layer 
of water which can enter or leave the boundary layer through a porous plate below. 

In Section 2 we calculate the exact solution for the basic state. In Section 3 we 
investigate the linear temporal stability of the flow to disturbances when the Reynolds 
number is finite. Since the basic flow is an exact solution of the Navier-Stokes equations, 
we are able to calculate the critical Reynolds numbers for a disturbance of arbitrary 
wavelength. By varying the viscosity and density ratios of the two fluids, we determine 
the stabilizing/destabilizing effect of the interfacial mode. We find that for both wall 
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blowing and suction, the interface significantly destabilizes the flow. More precisely, 
we show that the flow is susceptible to traveling wave disturbances at lower Reynolds 
numbers than is the case for flow in the absence of a water layer. 

The inviscid stability of a three-dimensional boundary layer was first comprehensively 
studied by Gregory, Stuart k Walker [9]. These authors used both experimental and 
theoretical techniques to develop an extensive understanding of the stability of the flow 
which forms on a rotating disk, and their findings have important consequences for the 
stability of general three-dimensional boundary layers. 

The experimental work of Gregory e/ al. [9] was based on the china-clay evapora- 
tion technique. They observed a regularly spaced pattern of equiangular spiral vortices 
which remain stationary, relative to the rotating disk. The angle made between these 
vortices and the radius vector of the disk was found to be in excellent agreement with 
the mviscid theory developed by Stuart. The prediction for the number of vortices was 
not, however, in such close agreement with the experimental observations. This dis- 
crepancy was attributed to viscous effects, and was resolved later when Hall [10] used 
a self-consistent asymptotic theory to study the problem. Hall extended the inviscid 
analysis of Gregory et al. taking into account non-parallel flow effects. His results were 
consistent with those obtained by the parallel flow numerical investigation of Malik [19], 
although this approximation is not valid at finite Reynolds numbers. In this work, Malik 
obtained a neutral curve for these stationary disturbances, and he also found a second 
stationary mode of instability which had been discovered experimentally by Federov 
Plavnik, Prokhorov k Zhukhovitskii [7], 

In Section 4 we consider the inviscid stationary modes of instability of the flow de- 
scribed in Section 2. Using numerical methods we calculate the eigenvalues and eigen- 
functions when the fluid properties are equated. We then compare these with our calcu- 
lations for air flow over a water layer. In addition, we use asymptotic techniques for the 
case when the fluids have similar viscosities. This gives a useful method for quantifying 
the onset of the interfacial instability. We find that stationary modes are susceptible to 
interfacial effects due to a discontinuity in the shear stress at the unperturbed interface 
position. In Section 5 we draw some conclusions. 


2 Formulation of the Basic State 

We consider the three-dimensional flow of two viscous, incompressible fluids above an 
infinite, horizontal, porous flat plate. The two fluids are immiscible and occupy separate 
regions. The upper fluid velocity is denoted by Uj and the lower fluid velocity by U^. 
We use cartesian coordinates, with the ( x*,z *) axes lying in a plane parallel to the 
plate which is positioned at a vertical height y* = —d. The porous plate allows us to 
model either the case of wall blowing, where there is a flux of fluid into the lower region, 
or wall suction where the normal velocity at the plate is in the -y* direction. The 
streamlines in the (x \y*) plane extend to infinity, the volume of fluid in each layer is 
then assumed to be constant and the interface between the two fluids is located at a 
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height y* = 6y* ( z*,t *) where if is an unknown function, and 6 is constant. 

The upper and lower fluids have viscosities /i x ,/i 2 and densities p x and p 2 respectively, 
so that the kinematic viscosities are v\ (= pi/pi ) and v 2 (= n 2 / p 2 ). We define the fluid 
velocity and pressure to be 


Pj=\ ,2 = Pj(y\z*,n, 

and the Navier-Stokes equations are 

dU* 2 dU* dU* 


— r U - T r • Q l f 

dt* 3 3 dy* 

dV* .dv: 


dz* 


dt* 

d J^l + vr- 

dt* 3 dy 


dy 
dW ; 


f- + "7 


dz* 

dW/ 

dz* 


dy* 



(1) 


(2) 

1 f 

x*pj dx* 3 

(3a) 

‘ dP* 

+ P jdy* ~ V > V V] ’ 

(3b) 


(3c) 

+ 2 - ° 

(3d) 


The Laplacian is defined as 


V 2 (-) 


d 2 (•) dMO 

dy* 2 dz* 2 ’ 


and the subscript denotes the upper and lower fluids respectively. The form of the 
velocity and pressure fields (1) and (2) corresponds to an exact solution of the momentum 
equations (3a-d), hence it is not necessary to make the boundary layer approximation 
when deriving the basic flow, and in the subsequent analysis. 

The tangential velocity of the lower fluid satisfies the no-slip boundary condition 
(x*U 2 , w*') (y* = = (0,0). The velocity perpendicular to the plate is prescribed by 

y* ^ y * _ _ ( {'j — y Q) where Vo > 0 corresponds to blowing, and Vo < 0 represents suction 

at the wall. 

The conditions far from the plate are given by 


IT 

u ' l ’ 


w; - Wo, 


as y —* oo. 


Define A = (pil/Uopi )* , l is a length-scale in the streamwize direction and the velocity- 
scale is VVo, so that 


(x *U\V\W m ) 


a (x,r,z), 

Wo(XU,V,W). 


II I 
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Time r, and pressure P/ =I 2 are made dimensionless by A/W 0 and pi 1T 0 2 respectively. 
We also define the following non-dimensional parameters 


Re 

m 


WpA pi 

/0 


/O 

/'i 


K 0 /?e 




1 / = — 
th ' 


7? e is the Reynolds number, k is the dimensionless normal velocity at the wall, and D is 
the depth of the lower fluid, scaled with respect to the length A. The parameters m, p 
and v are the viscosity, density and kinematic viscosity ratios respectively. 

At the interface between the upper and lower fluid layers, both velocity and tangential 
stress are continuous. The normal stress exhibits a discontinuous jump due to the effect 
of surface tension a. Using the notation 


[(•);]’= a -(•), , 

we obtain the following conditions which are applied at the non-dimensional interface 
position Y = Sfj 






= o, 


dw; 

dZ 


+ /o 1 


Wit = ml 

’s^) 2 ) (m+?n) f 

dZ) { dY^dZ 


dUj c drjdU ;] 2 

flj dY N dZdZ\ j 

2/q dVj 2 ^dr, (dW lM dV J \ 
V dZ) j R e dY R e dZ\dY + dZJ 


2/o dWj 
R e dZ 



— /«i J6 


d 2 fj 
dZ 2 



0, 

0, 

0, 

0, 


where J — o/p x W%A is the non-dimensional surface tension coefficient. In addition, we 
must satisfy the kinematic condition 


§ i ( y - w - 0 


-f(g + wM . at Y- SH. 


(4) 


Let us now regard the flow in each region as a small perturbation of the basic state so 
that with 8 < 1, 


(XV,V,W) M = ^ t L' W y {Y) + s (xu,V,W).(Y,Z,i) , 


(5) 
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and the pressure is written in the form 

i ftvrA y Pj - 

P,=u 2 l IWo ) + It 2 , ’’ 

Note that since S < 1, the unperturbed interface position is Y = 0. We substitute 
this flow into the Navier-Stokes equations (3a-d), and take the limit 6 -+ 0 to yield the 
following system of equations which determine the basic state. 


U x + V x = 0, 

vr+Otf-r.n-i = o, 
w; - v,w\ = o. 

U 2 V 2 — 0 ? 

= o. 

i>Wi-v,W 2 = o. 

The boundary and interface conditions become 


U 2 (-D) = W 2 (-D) = 0 , 

T?i(oo) = L, 

U 2 ( 0 ) = U\ ( 0 ) , 

V 2 (0) = 0, 

mV 2 (0) = U[ (0) , 


VA-D) 

Wi (oo) 
W 2 (0) 
Vi (0) 

mW' 2 (0) 


K, 

W, (0) , 
0, 

( 0 ) , 


(6a) 

(6b) 

(6c) 

(7a) 

(7b) 

(7c) 

(8) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


where (•)' denotes differentiation with respect to the normal coordinate Y. 

Before finding a solution to the above equations, we firstly analyze the behavior of 
the basic flow as Y -> oo. For large Y, the asymptotic form for V\ and W i can be 
expressed as 

X<1, 

L 2X 

2 


Vx = -( + r 0 x, 


W\ = 1 — A c 




exp 


(13a) 

(13b) 


1 3 15 , 

^ ~ £3 + £5 ^7 ^ 


1 + 


(-l) n (2n — 1) (2n — 3) • • • 3.1 

£2n+l 


+ 


(13c) 


where £ = Y + r , x — X (0 aiK l r » are constants ( see Rosenhead, Chapters V and 
VIII [26]). After substitution into equations (6a-7c) we integrate with respect to Y to 
obtain the following asymptotic form for x as Y °° 


X 


I [x" + ex'] , 


i + s + { 2 £ 


1 

2 


( e" 

£Sexpl-^- 


exp — 
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Having derived expressions (13a-c) above, we obtain numerical values for the basic flow 
using a fourth order Runge-Kutta scheme to integrate equations (6a- 7c) with respect to 
Y, from Y, x to -D, where Y^ is an arbitrarily large number. Initial values for t and T 0 
were chosen, and then improved in order to satisfy the no-slip conditions (8) at the wall, 
and the kinematic condition (11) to within a specified tolerance. For the case of a single 
fluid, (m = 1 = p) a step length of 1.0 x 10 -5 gave excellent agreement with the results 
published in Rosenhead [26] (chapter V on page 232). To model the flow of air over 
water we obtain a solution of the system governing the basic state using the viscosity 
and density ratios shown in Table 1 (see Batchelor [1]). The basic flow profiles U, V and 
IF are illustrated in Figures 1(a) and 1(b) for blowing and suction respectively. Each 
figure shows the velocity components with depth of water D = 0.5, 1,2,3 and 4. 

Given a constant depth of water D, we calculate the corresponding’blowing or suction 
K at the porous plate, the results are illustrated in Figure 2(a). For the case k > 0 (wall 
blowing), we see that the velocity at the plate increases almost linearly with depth of 
lower fluid D. With suction at the wall, the relationship is more involved. 

Before we examine the flow stability for both cases, we firstly discuss the basic flow 
properties. The way in which a constant depth of lower fluid is maintained may not be 
immediately obvious, especially in the case ofjvall suction. This is made clear by an 
analysis of the streamlines. By integrating XU x and XU 2 with respect to the normal 
coordinate Y we obtain Figures 2(b) and 2(c) which show the streamlines at a particular 
location along the spanwize direction. We have chosen representative examples: « = 
0.04, D = 2.0; and k = -0.12, D = 1.0. With a positive normal velocity at the plate, 
fluid enters the lower layer and moves towards a stagnation point at X = 0 = Y. The 
flow in the upper layer is directed towards the plate, in the —Y direction. This is a 
classical Heimenz stagnation point boundary layer solution together with an imposed 
crossflow W acting in the spanwize direction. 

For the case of suction, two stagnation points occur. There is a region above the 
interface (positioned at Y = 0), where V\ > 0 and If, < 0 as shown by the velocity 
profiles in Figures 1(a) and 1(b). The depth of this region increases with the suction 
so that with a depth D = 1.61 and /c = -0.1022 (suction is a minimum here) the 
two stagnation points almost coincide. For 0 > Y > -D the flow is towards the 
porous plate where the tangential velocity satisfies the no-slip condition. At Y = 0 
the kinematic condition (11) is imposed to prevent the transfer of fluid particles across 
the unperturbed interface (since the fluids are immiscible). The tangential velocity is 
continuous here although the gradient is discontinuous due to the viscosity ratio m^l 
(namely equation (12)). 

The relationship between k and D shown in Figure 2(a) can be analyzed as follows 
For k> 1 and Y ~ 0 (1) equations (6a- 7c) yield solutions 


U 2 

V 2 


Y + D 

pK 


+ 


K — 


Q' + P) 2 

2/9/c 
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W 2 = 1 - exp 


P(Y + D)k 

m 


Imposing the kinematic condition V 2 (Y — 0) — 0 yields 

D = (2 p)* k ~ 40.4k, 

which is in excellent agreement with the numerical values presented in Figure 2(a). 

For wall suction, the limit k -* -oo corresponds to the singularity in the depth D. 
A s K — ► — oo, we see that D < 1, and momentum conservation in the spanwize direction 
suggests the use of the scaled variable 

, (Y + D)\k\ P 


so that for £ ~ 0 (1) 


We now investigate the stability of the basic flow calculated above. We consider two 
distinct cases of physical interest: in Section 3 we look at the temporal stability of the 
flow when the Reynolds number is finite; in Section 4 we investigate inviscid stationary 
modes at high Reynolds numbers. 


3 Viscous Modes 

The aim of this work is to quantify the effect of the interfacial viscosity and density 
stratification upon the stability of the basic flow when viscous effects are included. For a 
single fluid, linear and nonlinear stability analyses have shown that unstable disturbances 
propagate along the attachment line. The three-dimensional basic flow is independent of 
the spanwize coordinate Z (and is therefore an entirely parallel flow). Hence we employ 
periodic boundary conditions (in that direction) on the flow disturbances. Such methods 
cannot of course be used for flows which are spatially growing (non-parallel) such as the 

Blasius boundary layer which forms on a flat plate. 

The flow described in the previous section is a first approximation to the boundary 
layer which forms on a swept wing, and is used to gain an understanding of the instability 
mechanisms which lead to transition from laminar to turbulent flow. To this end, we 
consider a convective instability in which disturbances propagate away from their source. 
For a discussion of absolute and convective instabilities the reader is referred to the 
review paper by Heurre & Monkewitz [13]. Following the work by Hall, Malik & Poll 
[12], we consider the temporal development of small amplitude perturbations having a 
normal mode expansion 

(ij,Pj=ia) = iv,Pj=h a)exp(ifc[Z-d]), (14a) 

(t7, V, W) 2 = (U, V, W) } exp (ik [Z - ct}) . (14b) 
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These perturbations are spatially periodic with wavelength 2tt jk and with speed c. 

The system of equations which govern the linearized stability problem are given 
by substituting equations (5 and Ua-b) into the Navier-Stokes equations (3a-d) and 

associated boundary/interface conditions and then discarding terms which are o(S) We 
obtain v ’ 


where 


(H 7 i) 

lh + V{ + ikWt 
L\{U 2 ) 

l\{v 2 ) 

L\{ W 2 ) 
U 2 4- V 2 + ik\V 2 


= 2U 1 U l + U' 1 V i +V,U[ 1 

= R'P[ + V l V l +V l V’, 

= ikR e P t + V X W[ + R'W\V U 
= 0, 


= 2 U 2 U 2 + U 2 V 2 + V 2 U^ 

R e P’ _ 

= -y- + V 2 v 2 + v 2 V', 


ikR e P 2 

P 


+ V 2 W' + R e W' 2 V 2l 


= 0, 


(15a) 

(15b) 

(15c) 

(15d) 

(15c) 

(15f) 

(!5g) 

(15h) 


L \(‘) = (')” “ ^’ 2 (•) — ikR e (iTj — c) (•) , 

Ll(-) = »{■)" -vk 2 {-)-ikR e (W 2 -c){-). 

The velocity perturbation to the lower fluid satisfies the no-slip condition at the 
plate, and the conditions at the interface are obtained by expanding the velocity and 
stress components as Taylor expansions about the unperturbed interface position Y = 0. 


m 


(*hv 2 , 


W 2 + t]W 2 = w x +rjW u 




Re 


+ W 2V + ikV 2 + W' 


2 ikr)V[ — „ 

+ W 1 r l + ikVi +W[, 


m 


t 77" > 

'£i + w 

R. + 2 


Re 

7/ 


P 2 + 


VPj 

PI 


m 


f 2 V V" . 2V' 


R] 


+ 


Re 


The kinematic condition (4) becomes 

7/ = 


| = 

Re 

= P 1 + 

— k 2 r)J. 

ReVi 


+ V[, 

rjP\ 2 V V; 

R\ Rl 


2v; 

Re 


ikR e {w,-c) -v[‘ 


(16) 
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The upper fluid velocity Ui must match the undisturbed flow as we move far away 
from the plate, we therefore require the perturbed flow to decay exponentially as Y 
becomes large. Ilall el al [12] showed that by replacing the basic flow by its asymptotic 
dependence for Y > 1 (equation (13a-c)), the perturbed velocity in the upper fluid has 
the form 

Ur ~ W\ ~ exp (-Y 2 / 2) , 1 M y _ QO _ 

Vi ~ cxp(-kY) , J 

Equations (15a-h) govern the stability of the lower fluid and are defined on the domain 
-D <Y< Si], whilst those for the upper fluid are defined for Si] < Y < oo. These 
equations in general require a numerical solution. 

For a given Reynolds number R e , and real wavenumber k, we obtain the correspond- 
ing complex eigenvalue c. The imaginary part, denoted c,, determines kc,, the linear 
temporal growth or decay of the perturbation to the basic state. When c; > 0 the flow 
is said to be linearly unstable and for c, < 0 it is linearly stable. 


3.1 Numerical Solution 

Solving the stability problem by means of a standard shooting method becomes pro- 
hibitively expensive as the Reynolds number increases. The rapidly varying nature of 
the eigenfunctions results in a loss of independence of orthogonal solutions due to the 
introduction of a ‘parasitic’ error at each integration step. High accuracy can only be 
guaranteed if the step-length is made vanishingly small. These difficulties were overcome 
by implementing a compact fourth order finite difference scheme of the form developed 
by Malik, Chuang k Hussaini [20]. This method was later used by Hall et al. [12] 
to investigate the attachment line stability of a single fluid, a detailed account of the 
implementation of this scheme is given by these authors. The method is applicable to a 
set of linear first order ordinary differential equations with an equal number of boundary 
conditions prescribed at each end of the domain. Our solution strategy is as follows, 
the equations describing the stability problem above have been formulated as two sixth 
order ordinary differential systems with coupled interface conditions. We define two 
column vectors 

T 

V’jsi.2 = (4>\ji 'foj, ^3j, <Sej) T = Vj, IT,, Pj,Uj, VFj) , 

where, as before, the subscript j- 1,2 denotes the upper and lower fluids respectively, 
and T denotes the transpose of the vector. The equations can now be formulated as 
twelve first order linear differential equations such that 


dY 

<P<hi 

dY 2 


6 


n— 1 

/ = 1,2, ... 6, j = 1,2, 

— ^2 i^ln) j <f>n j 7 

n=l 

1 = 1,2, ...6, j = 1,2, 
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bln = 


da 

Tv 


In 


"b 4 ^ (llpdpn- 


p=I 


Defining f 2 — k 2 mp 1 -f- ik(Wi — c)/? e , we find that the 6x6 matrix (<i/ n ) 2 has the 
following non-zero elements 


(" 15)2 = 

( fl 36) 2 = 

(043)2 = ikpV 2 R 

(« 5 i ) 2 = pm~ l (/ 2 + 2U2) , 

(062)2 = prn-'R'W^ 

(060)2 = pm~ l V 2 . 


?-i 


(021 ) 2 = - 1 , 
(« 4 i) 2 = pV 2 R~ e , 

(««)a = —771 7 ?£* , 
(052)2 = pm ~' U' 2l 
(oc 3) 2 = pm- 1 /- 2, 


(023)2 = — *&» 

(<>«), = -pR;' (/, + %) , 
(046)2 “ —ikmR- 1 , 

(055)2 = pm~'V 2, 

(064)2 = ikm~ l R e , 


The corresponding matrix (o/ n )j is obtained from the above by setting ?n and p to unity 
The numerical method is then derived using the Euler-Maclaunn formulae 

W = = <MK), 

= n-n-i, 


/?r 


= /»,/'#? #r‘\ ^ w , y 

J O l Vx”" ' I 1 n iTm * 


( 17 ) 

( 18 ) 


2 dT y 12 Ur 2 

+ 0 (/i*) . 

The nodes are distributed so that in the upper fluid 

Too + U 


dY 2 


( 19 ) 


9 1 = 

Y n = 


v ’ 

1 OO 

L x (n - 1) 


n = 1,2, ... + 1, 


iV<7i - (?r - 1)’ 

wheie N + 1 is the total number of nodes, is the edge of the boundary layer, and 
the scaling parameter chosen such that W x (Li/2) = 0 . 5 . Malik et al [ 20 ] showed that 
such a choice of Lx yielded maximum accuracy. Similarly in the lower fluid layer 

D + L 2 


92 = 

Y n = 


D ’ 

L2 (n - 1 ) 
Mg 2 - (n — 1)’ 


?! = 1 , 2 ,...,M+ 1 , 


such that W 2 (L 2 / 2) = k/ 2. 

For both the upper and lower fluids, equation ( 19 ) becomes 

*" - r-YE«4Ew-yE«;,-'r 

z P =i 1Z P = 1 1 p=i 


2 6 


- T? E ‘SrV;-’ = 0 . 

l ~ p=i 


n = 2 , . . . , N + 1, / = 1 , 2 ,. . .6, 
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which may be written in block-tridiagonal form so that the solution across each fluid 
layer is obtained efficiently. To this end, we introduce independent inhomogeneous ve- 
locity components at the interface, and equation (16) gives the corresponding interfacial 
deformation tj. We find a suitable linear combination of these three independent solu- 
tions, so that for a specified lower fluid depth D, Reynolds number R e , and wavenumber 
A', the conditions of stress continuity at Y = 0 are satisfied, and the complex eigenvalue 
c is obtained. When we equate the densities and viscosities of the two fluids, the numer- 
ical scheme yields exactly the eigenvalues found by Hall et al. [12]. When the imaginary 
part of the eigenvalue c is zero, there is no temporal growth or decay of the disturbance 
to the basic state, and the flow is neutrally stable. We then iterate to obtain neutral dis- 
turbances characterized by c, = 0. Figure 3 shows four neutral curves: an impermeable 
plate with k = 0; wall blowing with k = 0.137 and n = 0.4; and with suction k — -0.1. 
Inside the curves, c has a positive imaginary part and the perturbations {Uj,Vj,Wj,T]) 

grow exponentially in time. > . 

The eigenvectors given in Figures 4(a)-4(c) have been normalized so that the maxi- 

mum magnitude of each velocity component is unity. Figure 4(a) shows both real and 
imaginary parts of the three velocity components when the fluid viscosities and densities 
are equal. It has been verified that these (and other) eigenvectors are the same as those 
published by Hall et al. [12]. In Figures 4(b,c) we clearly see the discontinuities in the 
velocity and shear stress at the unperturbed interface position Y = 0 which is due to the 
difference in viscosities and densities of the air and water layers. It is this discontinuity 
which plays an important role in altering the stability of the flow. The neutral curves 
for air flow over water are drawn in the (k,R e ) plane. Figures 5(a) and 5(b) correspond 
to cases of wall blowing and suction respectively. These results are discussed in the 

following section. 


3.2 Discussion 

Before we discuss the novel results of our numerical calculations we firstly comment on 
the relevance of the exact solution, to the actual flow which forms on swept wings and 
swept cylinders. The boundary layer flow over a yawed, infinitely long cylinder was 
investigated and by Sears [27], (and in the unpublished work of Schubart). Their work 
is discussed in Chapter VIII of Rosenhead [26]. Using cartesian coordinates, the velocity 
components are expanded in powers of x/l, where z is the distance measured along the 
surface perpendicular to the cylinder generators and / is an appropriate length-scale. 
Close to the leading edge of the cylinder, and for a sufficiently large radius of curvature, 
the effects of curvature can be assumed negligible. The leading order solution (hig er 
powers of x/l are ignored) reduces identically to the flow which we have calculated in 
Section The accuracy of this approximation depends, therefore, on the geometry o 
the cylinder or swept wing. Results using this model will be most relevant to wing 

sections which have a flat nose. 

Since the basic flow is only a first approximation to the flow near the attachment 
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line, asymptotic methods based on a high Reynolds number assumption must be used 
to investigate the practical problem. 

It is worth making a few comments about the dimensional quantities in this problem. 
The velocity components in each fluid are made dimensionless using the spanwize free 
stream speed \Vq. The length scale A = is based on the streamwize velocity 

Uo and length /. In a practical situation then, the density and viscosity of the water and 
air would be fixed parameters (given in Table 1), as would the normal velocity at the 
surface, V 0 . We have shown in the previous section, that with a given value of k (the 
dimensionless parameter quantifying the amount of blowing or suction) we can calculate 
the corresponding nondimensional depth of water D. The actual height of the interface 
is therefore not a free parameter and is determined by the dimensional speeds U 0 , V 0 , W 0 
so that A is known and hence the depth d = D/A can be deduced. 

The results of our linear stability analysis are in excellent agreement with those of 
Hall et al. [12] when the fluid properties are matched across the interface (see Figures 
3 and 4). For a given wavenumber k, we calculate the Reynolds number which gives 
neutral stability. In the absence of suction or blowing, our numerical scheme yields the 
critical values ( R e ) c = 583.14, k c = 0.2881 in agreement with [12]. For R e < (R e ) c 
disturbances are damped and decay to zero exponentially in time. At points inside the 
neutral curve, the boundary layer is susceptible to traveling wave instabilities which 
propagate along the attachment line. 

An additional check on the numerical results is given by halving the step-size used 
in the finite difference calculations. Table 2 illustrates the accuracy of the scheme as the 
number of mesh points is doubled. 

For a single fluid (corresponding to the case when the fluid properties are matched), 
suction and blowing have opposite effects on the flow stability. As k(> 0) is increased, 
the critical Reynolds number decreases, and the flow is linearly destabilized by a smaller 
crossflow velocity. See for example, the neutral curves in Figure 3 with /c = 0.137 and 
K = and the results given by Hall et al. [12]. Suction however, can be a useful 
laminar flow control. The stabilization induced by negative normal velocity at the 
sui face increases the critical Reynolds number, as illustrated by the representative case 
K = '~0-l in Figure 3. We have also calculated neutral stability results for other values 
of k (namely k = —0.15, —0.2, —0.25). In each of these cases the flow is stable over the 
range 0 < R e < 1500 illustrated in Figure 3. The asymptotic results of Hall et al. show 
that as k — > — oo, ( R e ) c can be made arbitrarily large. This however, does not take into 
account the effects of nonlinearity. Hall k Malik [11] showed that solutions bifurcate 
subciitically from the upper branch of the neutral curve. The linearly stabilizing role of 
suction may therefore be destroyed by nonlinearity and transition may be enhanced by 
the unstable nonlinear modes. 

Upon introducing a viscosity and density difference across the interface, the results 
of the linear stability analysis are significantly altered. For the flow of air over water 
(the fluid properties are given in Table 1), we have obtained results in the case of both 
blowing and suction at the wall. With a positive normal velocity at the porous plate, 
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we have chosen the representative cases: k = 0.027, k = 0.04, and k = 0.137. These 
neutral curves are illustrated in Figure 5(a). To emphasize the interfacial efTect, we 
have also included the curve (broken line) corresponding to the neutral stability of a 
single fluid (see Figure 3). These eigenvalues were calculated by following the results 
given by fluids with matched physical properties, and gradually introducing viscosity 
and density stratification across the interface. As m and p increase, the interfacial mode 
destabilizes the flow. For any given Reynolds number, the band of unstable wavelengths 
is significantly increased. The upper and lower branches of the neutral curve open out 
and the critical Reynolds number decreases. For example, with k = 0.137 and D = 6.0, 
we obtain critical values k c = 0.499 and (/? e ) c = 97.81, whereas for matched fluids the 
critical values corresponding to k = 0.137 are k c = 0.309 and ( R e ) c = 315.12. 

With suction at the wall, the viscosity and density stratification across the inter- 
face also leads to destabilization and the flow is again unstable for a wider band of 
wavenumbers. In Figure 5(b) we show neutral curves for the cases k = —0.1026 and 
K = —0.208 which correspond to water depths D = 1.0 and D = 0.5 respectively 
(see Figure 2(a)). With a nondimensional water depth of 1.0, the flow is unstable 
even for small Reynolds numbers. Accurate numerical experiments yield critical val- 
ues ( R e ) e = 10.9887, k c = 0.7946. When the depth of the water layer is reduced (and 
consequently the suction parameter is increased) the critical Reynolds number increases 
along with the corresponding wavenumber. For example, with D = 0.5, we obtain 
(ft'} _ 82.0096, k c = 1.4410. It is clear then, that the usual stabilizing effect of suction 
at the plate, has been negated by the strongly destabilizing influence of the viscosity 
and density discontinuities at the interface. 

A comparison between the theoretical and experimental results is difficult. As dis- 
cussed earlier, in-flight calculations and wind tunnel experiments indicate that a water 
layer on the wing surface can have a detrimental effect on drag and lift. This is most 
likely due to the premature transition from laminar to turbulent flow. The interfacial 
traveling wave instability observed here is a possible contributing factor in this process. 
However, experimental investigations into the instability of superposed fluids have had 
limited success in quantifying the interfacial mode. Charles &; Lilleleht [3] and Kao & 
Park [18] studied the plane Poiseuille flow of oil and water in a channel. They found 
instability at large Reynolds numbers which appears to arise in the water layer and 
causes the interface to become wavy. It is not clear that this instability is caused by 
the interfacial effects, it is more likely that the presence of unstable Tollmien-Schlichting 
waves in the less viscous fluid (water) are being observed at the interface. This mode 
is present at high Reynolds numbers in the absence of a second fluid and is perhaps 
the one observed experimentally because it has the largest growth rate. More successful 
experimental results have been obtained for two fluid flows in cylindrical geometries, 
where traveling waves are often observed at the interface. The books by Joseph and 
Renardy [17] give a good review of recent experimental and theoretical investigations. 

The flow described here is a crude model of the actual flow of air over water on 
swept wings. To make qualitative comparisons between the theoretical calculations and 
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observable phenomenon would require a more sophisticated model in which nonlinear 
effects are taken into account. The methods adopted by Hall & Malik [11] could be 
applied to the two fluid problem in an analogous manner, although the nonlinear in- 
terfacial conditions would complicate the analysis. In addition, global methods could 
be used to calculate the complete set of eigenvalues, relating the interfacial effects with 
other modes of instability. 


4 Inviscid Stationary Modes 


We now investigate the stationary instability of two-phase flow of air above water over 
a swept wing, when the Reynolds number is large. As before, we regard the flow in each 
region as a small perturbation to the basic state. The normal coordinate must now be 
scaled on the Reynolds number, so that the fluid velocity and pressure are 


(XU, V, W) j=l 2 = (xUj, m Vj, Wj^j ( m) + 6 (xUj, V } ,W } ) (X, Y y Z ) , 
p 1 fU 0 XAR e 

Pj = h2 = ~ 2 V 


Wq 


2 ~n 


+P t +rp - 


( 20 ) 

(21) 


After substituting (20) and (21) into the non-dimensional Navier-Stokes equations and 
taking S -> 0, we recover the ordinary differential system (6a-7c) and boundary condi- 
tions (8-12) which determine the basic state. At next order the equations governing the 
linearized stability of the lower fluid layer are 


0(x»g,) afi 

2 dX +XV2 ^, 


1 I XR^V 9 ^ 7 i XW 9 ^ 2 i ^ 

gy +AH , 2— +- — 

Y It ~ L I !7 I f> 2 TT7 dV2 Re 81\ 

A + Vy + H-gp- + w 2JI + 

a + a vM + 


rTjW. 

dX 


dX 


dY 


dY 


dZ pdZ 

u 4 . 9 4 . ph ^ 

U2+ dX +Re dY + ~0Z 



), (22a) 


(22b) 


(22c) 

o, 

(22d) 


where 


V 2 (.) 


gj) ■ n d(-) , d 2 (-) 

dX 2 e dY 2 dZ 2 ‘ 


The corresponding equations for the upper layer are obtained by replacing p and v by 
unity in the above equations. 

Following the inviscid instability theory of Gregory ct al. [9] we expect the pertur- 
bations to the velocity, pressure and interface to have the following modal expansions, 
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with wavelengths scaled on the boundary layer thickness, 

(XUj, Vj, Wj) = {XV j, Vj, Wj) (Y) exp (^iRI J adX + PZ J , (23a) 

( 7 7 ,Pj) = (7/,Pj)(V / )exp ^iRi J adX + fiZj. (23b) 

In particular we consider a flow which is neutrally stable so that the wavenumbers a 
and are real. As R e -» oo, an inviscid zone will develop with depth 0 (R e y2 )- This 
inviscid region is asymptotically matched onto a viscous wall layer so that the no-slip 
conditions can be satisfied at Y = -D. By balancing inertial and viscous terms in 
equations (22a-d), we see that this viscous layer has thickness 0 yR e 2 3 J- The inviscid 
perturbations Uj,Vj, Wj and P, and wavenumbers a and ft are then expanded in powers 

of 0 ( R ; y 6 ) 

u } = u ]0 + r:«u j1 + ---, 

Vj = Vjo + Re 6 Vj\ + • • ' j 

Wj = w j0 + R;hvj] + -'-, 

Pj = Pjo + Re ® Pjl + • • • , 

_ L 

CV = Oo "t" Re 6 Ctj -}-■■■» 

p = p 0 + + -••■ 

Substitution of the neutral disturbances (23a-b) into equations (22a-d) yields the follow- 
ing leading order system of equations which govern the inviscid stability of perturbations 
to the upper and lower fluids when the Reynolds number is asymptotically large. 


iXVjUjo + XVjoU'j = (24a) 

tUjVjo = , ( 24b ) 

P 

ilfjWjo + VjoW'j = (24c) 

iooXUjo + Vj 0 + iPoWjo = 0, ( 24d ) 

aoXUjo + PoWjo = Uj o- ( 24e ) 


Eliminating U j0 , W j0 and P j0 from equations (24a-e) we see that Vj 0 
equation, (25a-f) in each layer. 


satisfies Rayleigh’s 


Vi (v " 0 - -ToUo) = CT'.V.o, 


ye [0, oo), (25a) 
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U* {K-^Vn) 

Ko(oo) 

Vto (-D) 

V 2O (0) 

K (0) 

Here Uj is the ‘equivalent’ two-dimensional velocity profile, and 7* = cvg + ( 3 $ is the 
‘effective’ wavenumber. Note that the continuity of stresses at the interface are satis- 
fied trivially in the limit as R e 00. The inviscid solution V20 is matched onto the 
viscous perturbation in the wall layer, and in view of the continuity equation (22d) this 
perturbation is 0 hence V 2 o satisfies the boundary condition ( 25 d). 

The point at which U j — 0 is denoted by Y = Voi and as Y approaches this value, [/ 0 

and Wq behave like 1 / ( Y — fo)- By careful choice of <Vo /Pa, U j is also made to vanish 
33 Y -* *o> so that V 0 has no such singularity, and a classical critical layer analysis is 
not necessary (see Hall [10]). 

4.1 Asymptotic Solution for Similar Fluids 

The above system may be solved numerically, to do this a suitable initial guess must be 
made for the eigenvalue 70. To assist the location of this eigenvalue, we firstly consider 
the analogous problem where the two fluids have equal densities, and the viscosity ratio 
is close to unity, that is m = 1 -f e, where t < 1. This case corresponds to the flow of 
two fluids with similar properties, this is a useful indication of the manner in which the 
interfacial effects can alter the stability of the flow. 

The basic flow and wavenumbers are then expanded in an asymptotic series as 

a oo + fOoi + • ■ ■ , 

Poo + e/?oi + • ■ • , 

Too + CToi + • • • , 

Q'ooO'oi + PooPox 
Too 

U jo + eU j\ + • • ■ , 

Vjo + eVji + • • • , 

W j0 + eW jl + ---, 

U jo + ft/ji + • • • , 

atooXUjo + ftooWjo, 

otooXUji + PoqW + qqxXU jo + AnH'jo. 


0-0 = 
Po = 
To = 

T01 = 

Uj = 
Vj = 
Wi = 

3 = 

Ujo = 
% = 


u 2 v 2 o, 

ye |-D,o], 

(25b) 

0, 


(25c) 

0, 


(25d) 

Ko (0) , 


(25e) 

K (0) + - 

— Hi) V,o (0) Z7^ 0 (0) 
mU 10 (0) 

(25f) 
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The leading order basic flow in each layer j — 1,2 satisfies 

y jo — v ]0 v j0 - (v j0 ) + 1 , 

w; 0 = 

V\o (°°) = 

I'K 10 (oo) = 1, 


V 20 {-D) = «, 

V 70 (-d) = o, 

W 20 (-D) = 0, 


with Vj,Vj, 7 ", Wj, and W], continuous at Y = 0. 
At 0 (e), 


Vn 

— 

F,o V", - 2V' l0 V' u + VuC, 

vl 

= 

V 2 oV 2l — 2V 2 o^2i + V 2 iV 20 — V- 

w" n 

= 

VioW n + f„w; 0 , 

W" 
VV 21 

— 

V 2 qW 2\ T VnW l0 — W 2 o> 

V'u (oo) 

= 

0 , 

Wu (oo) 

= 

0 , 

v n ( -D ) 

= 

0 , 

Wn(-D) 

— 

0 , 

V n ( 0 ) 

— 

0 = V 21 ( 0 ) , 

K ( 0 ) 

= 

( 0 ) . 

Vn ( 0 ) 

— 

K ( 0 ) + V ”o (0) 1 

W u ( 0 ) 

= 

W 21 ( 0 ) , 

W' n ( 0 ) 

= 

F' 21 (0) + vi/; o (0), 


To solve the above equations numerically, we require the asymptotic form of the 0 (e) 
correction to the basic flow as 7 -> oo. This is obtained in a manner similar to the 
derivation of equations (13a-c). We find that 


Vn = Ti + TiX - r x rox\ 

W u = - Q + S ) + r i A °) ex P • 

Where t u Ti and are constants to be found. The equations governing the basic state 
were then solved numerically and the results compared with the solution of equations 
(6a-7c), choosing a value of m close to unity. The results we obtained gave excellent 

agreement up to 0(e 2 ). 


Ill 
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The solution of (25a-f) may be obtained by solving the adjoint set of equations (see 
Coddington and Levinson [4]). We recollect that if M is an ordinary differential operator 
over a region N, the adjoint problem is defined by 

J n 'LA/ ($) dy = $M+ (<L) dy = 0. 

In our case the region N = [-D, oo), contains two sub-regions [-£>,0] and [0,oo). This 
however, does not present a difficulty, following the work of Blennerhassett [2] we define 
a vector 


Zj 0 < Y < oo, 
Z 2 - D < Y < 0, 


and a 2 x 2 real matrix S such that 


Sj 0 < Y < oo, 
S 2 - D < Y < 0. 


For upper and lower fluids (j — 1,2 respectively), Z j and Sj are then chosen such that 


Z; 


l * ) 


K 


Vifi 

Uj 



^7oo 


\ 


1 



Equations (25a-b) may then be written in vector form 7J = SZ, where V, 0 satisfies no- 
slip at the boundaries and Z is continuous across the interface. The adjoint problem is 
now defined by 


(Z + ) T [Z'-SZ \dY 


[(5 )r C D 

Z T [(z + )' + S r Z+J dY = 0 . 


Writing S + = — S T the adjoint system becomes 


(Z+)' = S+Z+, 
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where the adjoint function Z^ is also continuous across the interface, and satisfies 
no-slip at the boundaries. The problem is self-adjoint. 

We now perturb the viscosity ratio about m = 1, and write 


VjQ = Vjoo + cVjoi + • ■ ' i 

S_j = §jo T eSj, H , 

Z j = Zjo + cZji + • - ■ . 

Substitution into the Rayleigh equations in each fluid layer yields 


Zjo — 
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Vi 


J oo 


VjooUio 


V! 

Vi 


io 
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Vjoi 

Koo \V in Un - u ]0 u. 


^jioi ~" 


VinUio , Q°° ~ 

Vi 


j o 


=2 
Uj 0 


Sj 0 — 


I'jo 




JO 


7oo 


1/ 


jo 
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jo 


Sji = 


Vj\Ujo - U j0 Ui i 
=2 
Uj o 


2700710 
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UjoUji-UxU. 


jQ 


=2 

^0 


Neglecting. terms of 0(e 2 ), it follows that the momentum equations are 

0(1) : Z 0 SqZo 0, 

0(c) : Zj -SoZ, = SiZo. 


(26) 

(27) 
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Vectors Z 0 and Z x remain continuous across the interface and Vjoo and Vj 0 i satisfy the 
no-slip conditions at 1 — D and as Y — * oo. For equation (27) to have a solution, the 
forcing term on the right hand side, must be orthogonal to the adjoint function, hence 

/_ 0 D (z 2 + 0 )s 21 z 2 odv + ^(z+)s 11 z 10 dr = 0. 

After some manipulation we obtain 


2700701 ?o 
I 0 

h 

h 


y,io(o)V, 0 (o) 

=“ 1 \ — 1 2, 

U 10 ( 0 ) 

/ o roo 



V 2 
’ 100 

Ki'w-UuK 


=2 

Uu, 



I200 

U 2l U 20 -UnU 20 \ 

=2 

U 20 


dr, 

dY. 


(28a) 

(28b) 

(28c) 

(28d) 


The integrands of I\ and I 2 are regular since the singularity at Y — Yq is removable, 
for details see Coward [5]. We are now able to calculate 701, the 0(e) correction to the 
effective wavenumber, by finding a numerical solution to the leading order momentum 
equations (26) and the solvability condition (28a). 


4.2 Discussion 

The Rayleigh equations (25a-b) and associated boundary and interface conditions (25c- 
f) describe the inviscid stationary modes of the two phase flow with general viscosity 
ratio m and density ratio p. These equations were integrated using standard a finite 
difference method so that for given m and p the eigenvalue 70 was calculated to a high 
degree of accuracy. 

Figure 6 illustrates the dependence of 7 q upon the lower to upper fluid viscosity ratio 
for 0.8 <m < 24. The eigenvalue is a strictly increasing function for m > 0. The effect 
of density stratification is more subtle, since it does not appear explicitly in equations 
(25a-f), but manifests itself through the calculation of the basic flow. 

In the absence of a discontinuity in viscosity across the interface, the wavenumber of 
the inviscid stationary mode is 


7o = 7oo = 1-4899. 


Using the asymptotic methods for m — 1 <C 1, we obtain the leading order correction 
to 70 due to a small viscosity difference across the interface. The solvability condition 
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(28a) represents two simultaneous equations to determine unknowns o 0 i and /3 0 i (taking 
real and imaginary parts of (28a)). However, it is more useful to evaluate 

Ooo °01 + 0oo0o\ 


We find that 


To = Too + 2e7ooToi + 0 (c 2 ) , 

= 1/1899 + 0.1726e + O (e 2 ) • 


Figure 7 shows the value of yl evaluated using the numerical scheme for m close to unity. 
The broken line represents the calculation of Too + 2(m - 1)tooToi by the asymptotic 
methods described above. 

The eigenvectors illustrated in Figures 8(a) and^8(b) have been normalized so that 
their maximum values are 1.0. Figure 8(a) shows V w and V 20 when the two fluids are 
identical. We notice that the maximum velocity perturbation is at Y = 0.0839. Figure 
8(b) however, corresponds to the case m = 5. Although the velocity perturbation is still 
continuous, a discontinuity in the first derivative at the unperturbed interface position 
has developed due to the equation of tangential stress continuity. The maximum of Vio 
now occurs much further away from the inteiface, at 1 1.341. 

The orientation of the disturbances relative to the streamwize axis is determined by 


the wave angle $ such that 


Oo 

0o 


Ooo QqiAx) — &OO0OI 
0oo 0oo 

0.7514 + 37.88e + O (e 2 ) , 


+ 0 (/), 



For matched fluid properties the effective wavenumber and wave angle given above cor- 
respond to the single fluid case. As viscosity stratification is introduced, we obtain the 
above corrections to these quantities and these in turn are in agreement with our nu- 
merical results for general viscosity and density ratios. These calculations are based 
on an infinite Reynolds number assumption. This work could be extended to include 
viscous effects in an analogous manner to the method used by Hall [10] for the flow over 
a rotating disk. Viscous effects enter at 0 (/?7 1/16 ), the corresponding momentum equa- 
tions must then be solved to determine Uj \ , Vyi , Wj \ , . . . , and the solutions matched onto 
the inviscid flow. The analysis is, however, made more difficult due to the complicated 
interfacial conditions which match the flow across the two regions. 


5 Conclusions 

In Sections 3 and 4 we have considered both two and three-dimensional disturbances to 
the flow of air over water. The exact solution of the Navier-Stokes equations described 
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in Section 2, is a crude model of the flow near the leading edge of a swept wing during 
heavy rainfall. We have shown that the interfacial forces have a significant efTect on 
the stability of the attachment line flow. Viscous traveling waves are predicted at lower 
Reynolds numbers than is the case for air flow in the absence of a second fluid. The 
instability is due to the discontinuity in the viscosity and density across the interface 
between the two fluid regions and occurs with either blowing or suction at the plate. 

At infinitely large Reynolds numbers, the interface also alters the neutral stability of 
stationary modes of the form considered by Gregory et al. [9]. The three-dimensional 
basic flow is written in terms of an ‘equivalent’ two-dimensional velocity profile which has 
an inflection point when the velocity is zero. Consequently the critical layer is passive 
and the ensuing calculations of the eigenvalues and eigenvectors for three-dimensional 
disturbances follow in a straight forward manner. Using both general numerical methods 
and asymptotic techniques for the flow of similar fluids we have obtained the corrections 
to the disturbance wavenumber and orientation due to interfacial effects. 
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Figure captions 

Figure 1(a). Basic flow of air over water: velocity profiles with wall blowing and 
water depth D = 0.5, 1.2, 3, 4. 

Figure 1(b). Basic flow of air over water: velocity profiles with wall suction and 
water depth D = 0.5, 1.2, 3, 4.. 

Figure 2(a). Basic flow of air over water: water depth D and corresponding blow- 
ing/suction K. 

Figure 2(b). Basic flow of air over water: streamlines with wall blowing. 

Figure 2(c). Basic flow of air over water: streamlines with wall suction. 

Figure 3. Neutral curves with equal viscosities and densities: Impermeable plate 
K = 0, blowing with k = 0.137 and k = 0.4, and suction with k = -0.1. 

Figure 4(a). Neutral eigenfunctions U, F, W for a single fluid, R e = 119, k = 0.4. 

Figure 4(b). Real part of neutral eigenfunctions U, V, W for air flow over water, 
R e = 1580, k = 0.04. 

Figure 4(c). Imaginary part of neutral eigenfunctions U , V, W for air flow over water, 
R e = 1580, k = 0.04. 

Figure 5(a). Neutral curves: solid line corresponds to the flow of air over water with 
increasing wall blowing; dotted line shows the neutral curve for a single 
fluid with no wall blowing. 

Figure 5(b). Neutral curves: the stability of air flow over water with suction at the 
wall. 

Figure 6. Eigenvalues 70 as a function of viscosity ratio m. 

Figure 7. Eigenvalue 70 for similar fluids: a comparison of asymptotic and nu- 
merical results when the viscosity ratio m is close to unity. 

Figure 8(a). Eigenfunction: Equal densities and viscosities. 

Figure 8(b). Eigenfunction: Equal densities, viscosity ratio m = 5.0. 
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Tables 



Density 
g cm ~ 3 

Viscosity 

g cm~ l s~ l 

Kinematic Viscosity 

an 2 s~ ] 

Air 

1.225 x 10" 3 

1.776 x lO-'* 

1.450 x 10-' 

Water 

9.991 x 10" 1 

1.137 x 10- 2 

1.138 x 10~ 2 1 

Water/ Air 

8.156 x 10 2 

6.402 x 10 1 

7.848 x 10- 2 


Table 1 . Physical properties of air and water. 


N 

k 

${kc} 

K 

10 

3.300581 x 10 _I 

1.226919 x 10" 1 

0.0 

20 

3.378719 x 10" 1 

1.267951 x lO" 1 

0.0 

40 

. 

3.384238 x 10" 1 

1.270776 x 10- 1 

0.0 

80 

3.384613 x 10" 1 

1.270965 x lO" 1 

0.0 

160 

3.384638 x 10" 1 

1.270977 x lO" 1 

0.0 

10 

8.540221 x 10“ 2 

2.284554 x 10" 2 

0.4 

20 

8.428938 x 10" 2 

2.243776 x 10" 2 

0.4 

40 

8.415787 x 10~ 2 

2.238444 x 10~ 2 

0.4 

80 

8.414404 x lO" 2 

2.237832 x 10‘ 2 

0.4 

160 

8.414255 x 10" 2 

2.237762 x 10~ 2 

0.4 

320 

8.414239 x 10~ 2 

2.237753 x 10" 2 

0.4 


Table 2. Neutral eigenvalues with decreasing step-size: m = 1, p = 1 and R e = 800. 
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